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Abstract 

Consider non-intersecting Brownian motions on the line leaving 
from the origin and forced to two arbitrary points. Letting the num- 
ber of Brownian particles tend to infinity, and upon rescaling, there 
is a point of bifurcation, where the support of the density of parti- 
cles goes from one interval to two intervals. In this paper, we show 
that at that very point of bifurcation a cusp appears, near which the 
Brownian paths fluctuate like the Pearcey process. This is a univer- 
sality result within this class of problems. Tracy and Widom obtained 
such a result in the symmetric case, when the two target points are 
symmetric with regard to the origin. This asymmetry enabled us to 
improve considerably a result concerning the non-linear partial differ- 
ential equations governing the transition probabilities for the Pearcey 
process, obtained by Adler and van Moerbeke. 

1 Introduction 

Consider the probability that n non-intersecting (Dyson) Brownian motions 

Xt(t) < ... < x n (t) 

in K belong to a set E £ K, with all particles leaving from the origin at time 
t = and all forced to end up at b\ < b 2 < . . . < b p at time t = 1: 
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,...,b p ) 



( all 0^(0) = \ 

all Xj(t) G E for 1 < j < n n\ paths end up at b\ at t — 1 



y n p paths end up at b p at t = 1 j 

with Xli=i ni = n and with (local) transition probability 

p(t;x,y) := — ==e ~~ f~. (1.1) 
v 7rf 



A formula by Karlin-McGregor enables one to express this probability as an 
integral of a product of two determinants involving the transition probability 



(1.1) above. This further leads to a expression as (i) a GUE- matrix integral 
with an external potential, (ii) a determinant of a block moment matrix, with 
p blocks and (iii) a Fredholm determinant of a kernel. Finally it is also the 
solution of a PDE in the end-points of the interval E and the target points 
h, ...,b p . 

Throughout this paper, we shall be dealing with the case of two target 
points p = 2. In this paper, we show that, when n — > oo and when one looks 
through a microscope near a certain point of bifurcation, the non-intersecting 
Brownian motions tend to a new process, the Pearcey process, whatever be 
the location of the target points and whatever be the proportion of particles 
forced to those points. Tracy and Widom [28J showed this result in the 
symmetric case; namely when the target points are symmetric with respect 
to the origin and half of the particles go to either target point. Brezin and 
Hikami [9l [Hll [HI [12] first considered this kernel and Bleher-Kuijlaars jH] 
obtained strong asymptotics using Riemann-Hilbert techniques. 

The Pearcey process V(t) describes a cloud of Brownian particles, evolving 
in time according to a (matrix) Fredholm determinant, 

F v (all V(t 3 ) G E°, 1 < j < m) = det (i - (x K? X ) ) 

of the Pearcey kernel 

I(s < t) (x-v) 2 

' _g 2(t-s) 
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(1.2) 

The contour X is given by the ingoing rays from ±ooe i7r / 4 to and the 
outgoing rays from to ±ooe _J7r / 4 , i.e., X stands for the contour, all rays 
making an angle of it/ 4 with the horizontal axis. 







For s = t, the Pearcey kernel can also be written 



p(x)q"(y) -p'{x)q'{y) + p"{x)q{y) - tp(x)q(y) 



x-y 



;i.3) 



with 



q{y) ■-- 



2tt 



e 4-+ tJ f~ u vdU, p(x) :-- 



1 / i ' 



e 4 2 ^ v x 



dV, 



x 



satisfying both, the differential equations (using integration by parts) 

p"'{x) - tp'(x) + xp(x) = 0, q"\y) - tq'{y) - yq{y) = 0, (1.4) 
and the heat equations 



dt 2 P [ h dt 2 q{Vh 
whereas Kf t satisfies the following equation 

dKZ 1 



dt 



p'(x)q(y)+p(x)q'(y)). 



;i.5) 



1.6) 



The latter follows from taking d/dt of the kernel ( 1.2[ ), which has for effect to 



multiply the exponentials under the integral (j 1 . 2p with ^(U —V )/(U— V) 
\{U + V). 



Consider n non-intersecting Brownian motions, with < p < 1 and b < a: 



( 



a) 



p| {all xj (U) for 1 < j < n} 



Ki<m 



all Xj(0) = \ 
pn paths end up at a at t = 1 

[l—p)n paths end up at b at t = lj 



When n — > oo, the mean density of Brownian particles has its support on one 
interval for t ~ and on two intervals for t ~ 1, so that a bifurcation appears 
for some intermediate time to, where one interval splits into two intervals. At 
this point the boundary of the support of the mean density has a cusp. We 
show that near this cusp, the same Pearcey process appears, independently 
of the values of a, b and p, showing " universality" of the Pearcey process; see 
Figure [T} As it turns out, it is convenient to introduce the parametrization 
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p = - — 3 with < q < oo and let r := y' q 2 — q + 1. (1.7) 

Theorem 1.1 For n — > oo, the cloud of Brownian particles lie within a 
region, having a cusp at location (xoy/n,to), with 



(2. + a) 1 =l + 2 (r±-m\ (lg) 
q + 1 to \ q+1 J 

Moreover, the following probability tends to the probability for the Pearcey 
process: 

\l<i<»n ^ ^ ' ' J / 

= P 7 * ( D {P{n) n E = 0} 



v Ki<m 



;i.9) 



using the following constants 



; q 2 -q + l \ l/ \ - / tp(l-to) r(a-6) 
A* = > 0, , c := \/ =t ~~; > 0. 



g J V 2 g+1 

A = g 1/2 (a-x ) + g- 1 / 2 (&-x ) 
(a — b) 

In [1], Adler and van Moerbeke showed that the Pearcey transition prob- 
ability (1.11) satisfies a non-linear PDE, expressible as a Wronskian of the 
expression (1.12) with some partial. This was obtained from taking a scal- 
ing limit, when n — > oo, of the symmetric situation, i.e., where b = —a 
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and p = 1/2. It came as a surprise to us that considering the asymmet- 
ric case leads to a different non-linear PDE, when n —>■ oo, but nevertheless 



also expressible as a Wronskian of the same expression ( 1.12 ) with some other 



partial. A separate functional-theoretical argument then enables one to show 
that the expresssion (1.12) itself vanishes. This was one of the motivations 



1.1 



for finding the exact scaling as presented in Theorem 

To E = U[ =1 (?/2i-i, y2i) C K one associates two operators, a divergence 
and an Euler operator 



2r 

E 



2r 



d_ 

dyi 
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Figure 1: The Pearcey process for b = 0. 



Theorem 1.2 The log of the transition probability for the Pearcey process, 
which is non- stationary, 



,E) := logP 7 (V(t) HE = 0) (1.11) 

satisfies the following 3rd order non-linear PDE in t and the boundary points 
ofE, 



+ ^[^-24-2]dlQ-UdlQ,d E ^} =0, (1.12) 



dt 3 8 V dt / 2 B dt 



9 E 
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with "final condition", given by the Airy processn which is a stationary pro- 

— o /2 

cess, (by moving far out along the cusp x = 2 (|J ) 



iim ( m ^y { f' 2 n (-E) = 0) = det(/ - a; 



t^oo 



(-E) 



Remark: It is interesting to compare the Pearcey PDE with the Airy process 
PDE; namely for semi-infinite intervals E\ and E2, the 3rd order non-linear 
PDE for the Airy joint probability 

Q A (t] x, y) := logP- 4 (Aih) < *±^, A(t 2 ) < , for t = t 2 - h, 

reads 

d 3 Q A r 2 x d\ f d 2 Q A d 2 Q A \ f g'Q- 4 ^Q- 4 ) 

dtdxdy \ dx dy ) \ dx 2 dy 2 J \ dxdy ' ch/ 2 J 



1.13) 



with "final condition" 



lim (Afa) < Wi, ^(t 2 ) < « 2 ) = ^(mi)^^)- 

In the last section (section[7]), we develop -in a formal way- the central role 
played by the spectral curve (or Pastur equation [21]) in the steepest descent 
analysis used to prove the universal behavior of the kernel as iV — > 00 for 
the different problems of non-intersecting Brownian motions. The spectral 
curve is precisely the function which appears in the steepest descent analysis. 
The spectral curve associated to the problem provides the universal limiting 
kernel obtained after a proper rescaling of the variable around a singularity 
of the problem. 

1 The Airy process is a stationary process, which describe the statistical fluctuations 
of the process about the curve appearing in Figure [T] away from the edge and properly 
rescaled. Its probability given at any time by the Tracy- Widom distribution T(x). The 
latter is given by the Fredholm determinant det(I — A) of the Airy kernel A, restricted to 
the interval under consideration. 
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2 Non-intersecting Brownian motions on R, 
forced to several points 

In the expression below, 7i n (E) is the set of all Hermitian matrices with 
all eigenvalues in E. Note that in general one has the following, using the 
Karlin-McGregor formula^ (see (2D El EH QH H21 EEl E\): 



( allx j (0) = \ 

all Xj(t) G E for 1 < j < n rii paths end up at b\ at t = 1 



lim 

s a11 ™ -* ° j£« At 

Oil ■ • • . °n 1 — » °1 



111 ^ 



n p paths end up at 6 P at £ = 1 y 
det (p(t; Xj)) x<i j<n det(p(l-t; a*, <V))i<t',j'<n> 



» ni +^..+ rep _ 1 +i 
1 



, . . . , S n — t b„ 



n JH n tE 



\l *(!-*)) 



dMe-^ T ^ M2 ~ 2AtM UM 



x l +i e \ +blX dx 

E / 0<i<m-l, 0<j<n-l 

det 

\ \Je / Q<t<rip-1, 0<i<n-l 

det(J — H^) EC , (Fredholm determinant) 



/ 



(2.1) 



2 A„(xi, . . . , x n ) is the Vandermonde determinant. 
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where Hn\x,y) is the kernel (setting i fc = t e = t) 
Hj?\x,y)dy 



t k V , 2xV 
l-tu ' 1-tu 



— , dy fdvf du e 1 ^ n 



U-L 



tHl + Mf XX \V-br) U-V 

g l-t f l-i^ r= l 







for tt. > to 



for £ fc < ^ 



(2.2) 



where X is a contour consisting of the two incoming rays from ±ooe J7r//4 to 
and the two outgoing rays from to iooe - " 1 "/ 4 , provided no b r = 0. In the 
expression above, A t is the diagonal matrix 



A, :-- 



O 



O 



V 



I Tli 

t n 2 



with b,; = b, 



2t 



i-r 



I 



n.,. 



\ J 



The main expression appearing in ( |2.1[ ) contains the matrix integral 

' An(xi, ...,x n ) A 



(2.3) 



F n {E; bi,...,b p ) 



1=1 



n e 

n 



= _L /" (iMe -iTV(M 2 -2AM) ) (24) 
Ai ./«„(£) 

which is now viewed as a function of the boundary points of E and the target 
points bi, for which one assumes a linear dependence 

v P 
/Jcj&j = with ^^Cj = 1. 
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Introduce the following operators: 



d, 



(0 



sum of partials in the 
boundary points of E 

Euler operator in the 
boundary points of E 

p-i 



p-i 

5> 



d_ 
dbi 



d d 

Q ^ ~ ~ OIle Checks 9 b = ( 2 - 5 ) 



Proposition: [5] The expression logP n satisfies a non-linear PDE in the 
boundary points of the interval E and in the target points b{, given by the 
(near-Wronskian) determinant of a (p + 1) x (p + 1) matrix 



(Fx 



det 



\ 



f{ f' 2 f' z ... f; g 1 



F'l F» F£ ... F>> G 2 



\Fi 



(p) cb) z?(p) 



where the F e and are given b}|^] (G = 0) 



Ft 
Ge+i 

HP 



(df + Q<9 B )«9 £ logP n + n, 



H 



(i) 



o-(2)~ 

F j 



(2.6) 



(-c^e + ((e - 1)q + 2){df ) + Q<9j) logP„ + C t 
(l-e + 2b t d B ) (cf + c,d E ) logP n . (2.7) 



with C f = -2n £ 1(1 - c e )b e + Y, j7 t e b^b~ 
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Example: Setting He := < Hj> , Fg\ — < H\ ,Fi \ , one has: 

p=l det(p) ^=Ht = 

( F 1 F 2 
p = 2 =► F^det F{ F 2 f + f 



TP" T?" H i i H 2 
r l r 2 Fl F 2 



(2.8) 



= (H 1 F 2 + H 2 F l ){F l ,F 2 }' 
(H[F 2 + H' 2 F 1 ){F 1 ,F 2 } = 



3 Proof of Theorem 1.1 



Proof : It is easily computed by the Pastur-Marcenko method [21], which 
states that, given a diagonal matrix A = (ai, . . . , a„) and its spectral function 
dcr(A) := $^5(A — a^, the random Hermitian ensemble with probability 
defined by 

±- / dMe-^ M - A ^ } (3.1) 

%n Jn n {E) 

has, in the limit when n —* oo, a spectral density dz/(A), whose Stieltjes 
transform 

di/(A) 



/(*) = 

satisfies the integral equation 



X-z' 



3m z 7^ 0, 



d<r(A) 



A-*-t;2/(z)' (3 ' 2) 

and, in view of the spectral function da(X) := - ^ <5(A — a*), this becomes, 

n 1 / 



^ ai - z- v 2 f{z) 



Consider the situation (2.1), where n, particles are forced to fej. Then, upon 
setting the variance v 2 = 1 and defining g(z) := f(z) + z, the equation (3.3) 



rea> 



9-z 



i 9 ~ h 



0, with £j 



71; 



n 



(3.4) 



See Section 7, formula (7.18), for further comments on the "fraction numbers" £j 
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The density of the equilibrium distribution is then given by 



dv{z) 
dz 



1 



7T 



f(z)\ 



1 



7T 



3m g(z) 



with z G 



Remark: It is precisely this Pastur-Marchenko equation (3.4), which will 
appear in the argument of steepest descent for the corresponding kernel in 
section 7. Namely, equation (3.4) is the derivative (7.19) of the S-function, 
appearing in (7.17). 

We will now specialize to two target points p = 2, where ni = pn and 
ri2 = (1 — p)n. From (2.1), it follows that for non-intersecting Brownian 
motions forced to b\fn < a\/n, one has 



e -%Tr(M-A t r dMj ( 3 _ 5) 



n JH n { E. 



t(i-t), 



with 



A, :-- 



a 



\ 



O 



a 



1 ' l-t 



O 



V 



pn 



1(1 



(3.6) 

For A 4 as above, the integral equation (3.2) becomes the algebraic equation 

P 



1 — p 

g-z+ — ^+ 



o 



9-13 g-a 
which, upon clearing, leads to a cubic equation for g, 

G{g) := g 3 -(z + a + (5) g 2 +(z(a + (3) + a/3 + 1) g-(a/3z+(l-p)a+p/3) = 0, 

, , (3-7) 

with roots given by g = q- 
inant in z, 



/ Ai + y/r — with a quartic discrim- 



Z - Zi 



0; 
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q, r are polynomials of z,a,(3. Thus one finds the following 

|3fm g{z)\ 



| I _ \^sm g(z)\ for z such that Ai(z) < 



dz ix 

for z such that Ax (V) > 

Therefore the support of the equilibrium measure will be given by 

either two intervals [z 2 , z ] U [21,23] 

or two intervals touching [z 2 , Zq] U [zo, Z3] 

or one interval [z 2 , Z3) 

for the real roots of Ai(z) = (a — (3) 2 Ylo( z ~~ z i) = 0- Thus depending on the 
values of the parameters a, (3 and p, there will be four real roots or two real 
roots, with a critical situation where two of the four real ones collide, say 
Z\ = Zq- The latter occurs exactly when the discriminant A 2 (with regard to 
z) of Ai vanishes, namely when 

A 2 (a,f3,p) = 4p(l-p)p((p - l) 3 - 27p(l -p) P Y = 0. 

V / p=(a— /3) 2 

This polynomial has one positive root (the others being imaginary), and one 



checks that, taking into account a > (3, and defining p and r as in (1.7) 

1/2 q+l 



a: 



(3 = ^f-p= (3P 1 / 3 (1-pf 3 + 3p 2 / 3 (l-j9) 1 / 3 + l) = ^ > (3.8) 



For this precise value of the parameter a— (5, two of the four boundary points 
of the support coincide, namely the roots zq and Z\ of Ai(z) coincide: 

n 2q - 1 q-2 , , 

z = Zl =f3+— = a + - . 3.9 

r r 

This double root z is found by stating that, under the condition a — (3 = 
(q + l)/r, the polynomials Ai and A[ have a common root; in other terms for 
appropriate choices of ci, 02,03, some linear combination of the polynomials 
Ai and A[ becomes a linear polynomial (4z+ci)Ai(,z) — (z 2 +c 2 z+C3)A[(z) = 
c^z + z 5 for some C4 and C5. Since z is a root of the left hand side, it also 
must be a root of the right hand side. This is to say z = —c^/c^, yielding 



the expression (3.9). 
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The critical time to is then obtained from setting t = to in (3.6), 



a = a 



2t ° and P = bJ 2t ° 



1-t 



1 - to 



(3.10) 



from which one computes t , by taking the difference and by using (3.8), 

, = 0?+l) 2 

(g + l) 2 + 2(a-o) 2 (g 2 -g + l)' 

and from which one further computes 

t I 2t q + 1 



c 



*o(l -*o) , r(a-b) 

o = *o ~r and — 

2 q + l c 



1 — t ( a ~ b)r ) 



(3.11) 



confirming the expression for cq in (1.10). Then from (3.10) and (3.11) one 
deduces 

to a{q+l) _ Jo b(q+l) 



a = a 



Co (a — b)r Co (a — b)r 



Defining x := z c , one computes from (3.9), 

x (2a - b)q + (26 - a) 



Co 



(a — 6)r 



(3.12) 



(3.13) 



Next, one computes the double root of the G-equation (3.7) for the value 
z = z . Indeed, using (3.9) and (3.12), one checks that for some root g = g , 



dg dg 2 



where 



^ := 3 {Z0 + a + (3) = 3 i2a 7 +b 7 + — ) = T^Wr 



and from (3.12) and b < a that 



/? < g < a. 



(3.14) 



(3.15) 



The point Zq in (3.13) refers to the matrix integral variables on the right 
hand side of (3.5), which can then be transformed into the Brownian motion 
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variables, according to E — E^ t<yl 2 ^ at t = t , which gives the transfor- 
mation from the matrix integral variables to the Brownian motion variables. 



and (3.13)) 



Hence the critical point in the Brownian picture takes place at (using (3.11) 

(2a - b)q + (2b - a 



{x y/n,t ) = (z c y/n,t ) 



q + 1 



inVn, tr 



To prove the main statement (1.9) in Theorem |1.1| start with the kernel 



(2.2) of the non-intersecting Brownian motion: 

H n (x,y;t k ,U)dy 

dy 



2ttV(1-4)(1-^) 



dV 



tV 1 j 2xV 



dU 



L+iU 



W 1 | 2yU 



U -b 
V-b 



»2 



U -a 



»i 



V-a U-V 



One first needs to prove that for some ip n (\, r): 
lim ip n (x, t k )H n (x, y; t k , U)<p n (y, 

n^oo 

= K r (^r];T k ,T(,) 



rcscaling 



4vr 2 



du v 



x 



du u 



e 4 



e 4 



with the rescaling j x y j 



U = t + (c fi)' 



2t, 



n 



1/2 : 



c z n 



1/2+ A 




with constants A, /i, given by (1.10). 



(3.16) 



(3.17) 



Consider the change of variables U := c °" v/ " . The form of the Brownian 

to 



motion kernel (2.2) suggests, by putting the two ^-factors of the integrand 



in the exponential, the function F(u), which one observes, at leading order, 
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is closely related to the function G(u), defined in (3.7); nameljj^] 



F(u) 
with 



uz + p\og{u — a) + (1 — p) log(w — (3) 



•I! P C| ) 



F'(u) 



G{u) 



(u-a)(u-f3) a= ato 



Remember from (3.14) that 



uo := go 



aq + b 



(a — b)r 

is a root of G(u) = and two of its derivatives, 

G{uq) = G'(uo) = G"{uq) = and G"'{uq) = 6 
and then one computes, since^] (uq — o){uq — (3) = —q{q 2 — q+l] 



(3.19) 



-i 



F'(u ) = F"(u ) = F'"(u ) = and ^F^(u ) = - ^ 



(3.20) 



and so 

F(u) = F(u )-^(u-u ) 4 + O(u-u ) 5 . 

In the calculation below, the first equality = is obtained by doing all the 
substitions below, except for the last one u i— > u u , whereas the second equality 
= is obtained by the substitution to the new integration variable !ihu u ; the 
expression after = contains a term n 1 / 4 ^, which contains the new integration 
variable and which blows up as n 1 / 4 . Hence this coefficient must be put = 0, 



5 Thc function F(u) should contain the term log v ^ c ° ; however as the same rescaling 
is made in the v- variables, this same term will appear with a different sign and therefore 
they will cancel. Consequently, this term will be omitted. 



upon using (3.9 1, a — (3 = (q + l)/r and the root u = h(zo + a + (3) as in (3.14) 
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yielding the value of A as in (1.10). In the next equality one uses this value, 
thus yielding in the end, 



tU 2 2Uy 



+ n 2 \og(U-b) + ni log(U-a) 



rii = np, n 2 = (1 — v) n 

jj _ c u^/n 



to 



t = t + (c ^ 2 ^ 



a i — > a^fri, b t— > b^/n 

y = c (zon^+ATE+n^) 



tiF{u) + n 1//2 TgU\^ 



u to A 

x 

2 c \i 2 



, , 4 2 / M *0 ^ 

+t n/i r € x 

V 2 c [i 2 



+ 0(n- 1 / 4 ) 



tiF(uq) + n 1 / 2 tiUqh 2 



u t -4 

7; n 

2 c /i 2 



+n 1/4 (r^ u (M - — x ^) - r/M )/i + ^K 1 ) 

c /r 



+0(n- x / 4 ), 



1/4 

72 7 UofiT] — 



using the value (1.10) of 



to j,, \ g 1/2 (a-x ) + q~ 1/2 (b-x ) 



A=y/q+^ (b-x ) 

Co 



(a -b) 



Similarly, 

tV 2 2Vx 



1-t 1-t 



+ n 2 logtV-fr) + rii \og(V-a) 



ni = np, ri2 = (1 — p)n 



t = to + (co/i) 2 J^ 
a i — > a-y/n, 6 i — > b^fn 

x = c (z ™ 1/2 + ^Tfc + At^i) 



f = v + 



1/4 
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2, ,2 



+0(n- 1 / 4 ) 



Moreover, using 



2d 



— 1, together with the rescalings above, one finds 
dU dV drj dw u duj v 



2n 2 y/(l-t k )(l-t e )U -V 4vr 2 (u u - u v ) 
Summarizing, one obtains: 



+ 0(n- 1 / 2 ) 



-dy 



27T 2 y/(l-t k )(l-t e ) Jc J L+ iR 



dV 



dU 



_*JsYl + 2xV_ 
C 1 ^ t k 1 ~ t k 

t t U 2 . 2yU 

e 1 ~ t e 1 ~ t t 

U-b\ ri2 fU-a^ ni 

x 



V-b 



dw„e~ a r+ 



4 2 



1 



+ 0(n- 1 / 4 ), 



where we deformed the u, v contours, by translating them by Uq, so they are 
no longer emanating from 0, but from u - Moreover, taking into account the 



extra-piece appearing in (2.2) for r k < ti, one computes 



dy 



y/n(tt - t k ) 



e t i~ t k e 1 ~ t k 1 ^ t t 



U = t + (c n) 2 i&2 

x = c (zon 1 / 2 + Ar k + (i^z) 

y = c {z a n x l 2 + At £ + vjji) 



dr\ 



(i-v) 



^F^T (l + 0(n- 1/4 )) . 



y/2n(T t - T k ) 

In other terms, setting 

D(£, t) := diag (. . . , e -«oA** 1/4 e -iV^^ a e -^to«g^ j _ _ 



(3.21) 



Kk<m 



19 



one has the following limit, upon using the (£j, x, y)-rescaling in (3.21) and 
upon conjugation of the matrix kernel, 



hm D(^T)(K ik fy(x 1 y)) 1 ^< m D(ri,T)- 1 = {K*„(t,r,)) 



which leads to the desired kernel (1.2), upon replacing the integration vari- 
ables uj u — > U, io v — > V. 

Since the above argument is obviously formal, one needs to make a rigor- 
ous steepest descent analysis on the conjugated kernel above. In the follow- 
ing section, steepest descent contours will be found, depending on whether 
q < 1, q = 1 q > 1, or what is the same p > |, p = ^, p < |. Notice the 
duality q < 1 q > 1. 



4 Steepest descent analysis 

In this section, we will deform the contours for both the u and v integration 
into steepest decent contours. They will be as depicted in the picture below; 
all lines are at an angle of 0, it /A or 7r/2 with the horizontal line. 

After having set U :— CoV ^ 1 and V := CQV ^ ri in the exponential appearing 
in the kernel, one was led to a function 



F(u] 



u 



uz +p\og(u — a) + (1—p) log(w - /?), 



where (upon using (3.9), (3.12) and (3.19)), 



2q-l ,9-2 at bt q 

z Q = — = p H = a H , a = — , p = — , u = p + -. 

c r r c c r 



Then one checks (with p = (1 + q 3 ) 1 and r = v/g 2 — q + 1) 



Uq — a = — - 



1 



r 



u - (3 



(4.1) 
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q = l 



q>l 



q<l 




(i) First we derive the steepest descent contour for the g-independent u- 
integration, which is the vertical line through u . Indeed, one checks that 
for 

MeF(u + iy) = ^(u 2 -y 2 )-u z +^ \og((u -a) 2 +y 2 ) + (l-p) log((u - (3f '+y 2 ) 
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and using (4.1), the derivative equals 




dy 



$teF(u + iy) 



y 2 + ^ 



3 y i r 2 _ 3 j a positive 

~" (y 2 + ^)(y 2 + £) ~ _2/ ' function 



showing KeF(M + has a maximum at y = 0, which takes care of the 
M-contour. 

(ii) The next point is to deal with the ^-integration. Returning to the 



non-intersecting Brownian motion kernel (|3.16|), the F-function goes with a 
sign on the f-contoi 

$teF(u + x(e + i)) 



negative sign on the v -contour and thus (with e = ±1) 

1 



((u + xef - x 2 ) + (u + xe)z 
log((w — a + xe) 2 + x 2 ) 
log((w - (3 + xe) 2 + x 2 ) 



2 
P 
2 

1 — p 



u — s(l — i) — X 




Uq — s(l + l) — X 



u + s(l + i) + x 



u Q + s(l — i) + x 



0<q<l 



1 < q 



r\q — 1| 
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Using (4.1),one checks for e = ±1 



^(w — a + xe) 2 + x 2 j (Juq — (3 + xe) 2 + x 2 ) -^$teF(u + x(e + z)) 



for 



4x 3 



y - (q - 1) ex + 



_0 
<9x 
Q 
r 



—x x 



-oo < X < 



r(g-l)' 



a positive 
function 



e = 1 



r(g-l) 



< X < OO, £ = — 1 



-OO < X < OO, £ = 1 

-oo < x < co, e = — 1 



if g > 1, 



if q 



(4.2) 



r(l-g) 
OO < X < 



< X < OO, 

5 



1 



-1 



if q< 1, 



r(l-q)' 

Therefore the function — 3fJeF(M + x(e + i)), restricted to the segments spec- 



ified by (4.2), has its maximum at uq. One then completes those segments 



by horizontal lines starting from the end of those segments as in the figure 
above. Along those horizontal half lines, one must check that the maximum 
is attained at the points uq ± s(l — i) and uq ± s(l + i). To carry out this 
computation, the four horizontal segments can readily be represented by 



uq + 6s(e + i) + 5ex 



with 



5 = -e = l 




-8 = e = 1 / x - 5 = - e = l 

in the four corresponding regions of the figure above. 
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To deal with the horizontal segment, setting s = , 9 T , , one finds 

° ' ° r-| g — 1 1 ' 



-9feFO + Ss{e + i) + Sex) 



--((u + eS{s + x)) 2 

_P 
2 

1 — p 



s 2 ) + (u + eS(s + x))z 



log ((w - a + e<5(s + x)) 2 + s 2 ) 

log ((m -P + e5(s + x)) 2 + s 2 ) . 



One computes, 

(uq + (s + x)eS - a) 2 + s 2 ) ((u + (s + x)e5 - (3) 2 + s' 



d_ 

dx 



F(u + 8s(e + i) + Sex) 



'r(q-l) 



r 5 (q — 1) 



qz 5 + eS(q 2 + 3q+ l)z 4 + 3(q + l) 2 z 3 + eS(q + 3)z 2 
+4(g 2 + q + l)z + eS(l + q 2 ) 



.r(q-l). 



(4.3) 



Since x > is increasing as one moves away from uq along the horizontal 
lines; one has simultaneously, 

--—2— eS = +l, q-l>0, z>0 



r(g-l)' 
g 

r(g-l)' 



eS = -1, q - 1 < 0, z<0 



and thus, since g > 0, the right hand side of (4.3) is > 0. Thus the above 
derivative is negative on the four lines and so, when one moves away from 
uq along the horizontal paths in all four directions (as in the picture above) 
the function — 9fte F(z) goes down so that combining both calculations, the 
maximum will be attained at u . 

In order to show that in the limit the Pearcey kernel is obtained, one 



along the contour in a neighborhood of radius n 



-1/4 
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Then 











? 


ft 



< 



S < nao ~ £ , as n —>■ oo. In this range lemma 4.1 will apply, whereas outside 
this neighborhood, the rest of the contour makes no contribution, because of 
the steepest descent estimates. So, one needs the following estimate: 



Lemma 4.1 Given the function F(u) as in (3.18), one has the following 
estimate, 



n 



< 



645 5 
5nV4 



q + 
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Proof: By Taylor's Theorem and using F'(uq) = F"(uq) = F'"{uq) 
and \F^ v \uq) = -g, as in (|3.20|), one has 



1 f 5 



n 



4! Vn 1 ^ 



< 



n 



1/4 



max 



From the explicit expression (3.18) for F, from the fact that Uq — f3 = q/r 
and a — uq = l/r and that (3 < u < a, one deduce^ 



max 

\u-u \<- 



1/4 



U 



5! 



< 



< 



< 



1 



sup 

l«-«0|<^74 



1 — P 



+ 



P 



(w-/3) 5 (w-a) 5 
2 



5min(|a - « - ^ | , |«o - P - ^ji 



Lh5 



64 



5 \(min(l,g)) 



by picking n large enough 
such that ^74 < ^ min(l, g) 
since 5 < n^^ 6 



< 



64 



q + 



ending the proof of Lemma 4.1 



5 Proof of Theorem 1.2 



In this section, we denote by a and b the Brownian motions target points, 
where we put 6 = 0. We denote the old time in the Brownian motion formula 
by t and the new rescaled time and space in the Brownian motion formula 



(Theorem 1.1) by f and fj. Set E c = U[ =1 (?/ 2i _i, y 2 i) C E. Let x be the 
spatial variable for the matrix integral and a the variable appearing in the 
diagonal matrix, the other one being = 0. 



7 Remember r = y/q 2 — q + 1 
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The reader is reminded of the different players in the argument below, in 



accordance with (2.1) 
F ST (t,y,ay/n) 



all Xj{i) G£for 1 < j < 



n 



all 2^(0) = 

rii paths end up at ay/n at i — 1 
ri2 paths end up at at i = 1 



and, setting b\ = a, b 2 = 
1 



F n (a,Xi):= — / A n (x u ...,x n )Y[A ne (x {e) )Yl 



with P n (a,Xi) satisfying the PDEQ as in (|2J( and (Q 

/ F 1 F 2 \ 



det 



V 



0. 



with Fi and ifj given in (2.7). From (2.1), one has the following relationship 



F Br (t, y, ay/n) = P n ay/n 



2t 



l-V Vl V t(l - 1) 



FJa,Xi 



and also from Theorem 1.1 



F Br (t,y, ay/n 



t=t +(c fi) 2 , y=x n 1 / 2 +c AT+c ii-jlj I 



= f v (v{f) n UL^fhi-ufhi) = 0) + o^- 1 / 4 ) (5.1) 

Setting 6 = 0, the formulae (1.7), ( 1.8[ ) and (1.10) in Theorem 1 1 . 1 1 simplify: 



2t q + 1 2q-l /t (l - t Q ) at r 

Xq — — — —at , Co 



1 — to ar 



q + 1 



q+1 



q-{2q- l)t y/F 
A = , n 



Q 



1/4- 



(5.2) 



8 Given E c — U^ =1 (x2i-i 5 x^i) C M, the prime in the formula below denotes ' := ° 



dxi 
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For convenience, we set 



r 7] 
r := — , r}:=— z := 

V 7 ? 9 1/4 



g 2 - g + 1 



n 



1/4 



One finds, using z 2 = r/^/n and formulae (5.2) 



t = t + (c aO : 



2f 



n 



*o 1 + (1 " *c 



r f q 2 - q + 1 



to (1 + (1 - t )rz 2 ) 



1/2^ 



Moreover, one computes 



a = a-wn 



2t 
1-i 



na 



2t l + (l-t )rz 2 



1 — to V 1 — toTZ 2 



q+1 l + (l-t )rz< 



1 - t TZ 2 



We also have, using the formulae (5.2) above, 

77 



t(l-t) C ° Vc 



;r 



n + At + /i- 



n 



1/4 



t(l-t) 



c (2? - 1)^- + (q- (2q - l)t )^ + 



t r] A/r 



v/a + a-to)^ 2 )^!-^ 2 ) v z 

One also checks 



q l/2 1 g l/4 n l/4 



2g-i , 

— + (q-t {2q-l))r + V z 



n 



n 



n\ = np 



n,2 = n(l — p) 



q 3 z 4 
(9+1)' 



g 3 +l (g+l)(g 2 -g + l) {q + 1) 
Consider the map T z : (r,r]i) i— > 
:= T z (r,r]i) 

\ q+ 1)^/1 + (l~t )rz 2 { 2 -^ + (q-(2q-l)t )r + m z) 
zV(l- *orz 2 ) ' + (1 - to)rz 2 ) (1 - t T2 2 ) 
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and its inverse T, 1 • 



a 2 z 4 - (q + lf 



a(xiq - aq + Xi)z 4 - (g-l)(g+l) 2 

^PWifH))' z 3 (t z*a 2 + (q + i) 2 (i-t )) 



Then setting 



logP n (a,Xi) = F(r,T ]i ) = F(T- 1 (a,x i )), 



setting B = UJL-^a^i-i, %2i) C R, taking the derivatives 9 a := ^2 ^f:> e ^ := 
^XjT^- and then taking a series in z, the functions Fj and ifj in (2.7) have 
the following form, 



= ~9.bgP„ + ni = ^— 
da q + 1 



: -4 + z -2 F //_ 22 -i^! + 0(1) 



d 1 

— + 9 )d logP n + n 2 = 

'9a xJ * 6 g+1 



g 3 z -4 + qz~' z F" + 2z 



Or 

BF' 

It?" , o~-l _ 

<9r 



+ 0(1) 



* = {^,F 1 } a + {^),F 1 } ( 



(<?+l) 4 
2q 3 



q{2q 2 + 3q + 2)F'"z- g 
-{4q 3 -l + 6q 2 + 3q)^fz^ 
+0(z- 7 ) 

-q(2q 2 + 3q + 2)F"' 'z~ 9 
+ (Aq 3 -l+6q 2 + 3q)^fz~ s 
+0 (z- 7 ) 



with 



(2) 



H. 



(2) 



(i) 



(i) 



(l- e-+ aA + 2aa.)(-A)togP B 
(l-s. + aA)(A + a.)l«gP B 

-2^-logP n -2m(a+— ) 
aa a 

1 + e. - a— — logP n + 2^ 
aa aa a 



Then one computes (2.8), setting ' = d x and setting 8e = ^2 ^ for E 
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\ 


Hi 


+ 






Fi 


F 2 J 




H[ 


+ 








F 2 J 


J 



(a,Xi)=%(T,r)i), With 



V9' 



Vi 



-74 



6+1/2 



q-1 



(9 + 1) 



5 E 



9fi 



-18 



+ 



8 \ \dt 



d 



dl\og¥ v ^\ +0(q- 1) I z- u + 0(z~ lii ) 



where, setting Q(t, E) : = logP p (V(t) HE — 0), 



X := 8 



e,,-2l d - 2 I - ' N/rO.^ 



<9t 



where we made use of (5.1 ), which states that logP„(a, Xj) = Q(t, E) + 0(z). 
• For q 1, the function logP^, which is independent of q by the universality 
result, satisfies the differential equation, given by the leading term z~ 18 , 



{<9|logP p ,X} =0. 



(5.3) 



• For q = 1, the z 18 -term vanishes and thus logP satisfies another equation, 
namely the one appearing in the z~ 17 -term, 



al|io g p^,x} = o. 



(5.4) 



This means that logP 75 satisfies the two equations (5.3) and (5.4). 

Thus for E = (x, y) C R, setting u± = |(y ± x), and logP 7 "' = H(t; \(y + 



<9 B logP p 
e E \ogF v 



ox ay ou + 

(x— +y— )\ogF v 
ox oy 



9 9 ^ u 

[u+- h U-- — )H. 



du. 
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Since the wronskian of two functions equals the derivative of the ratio, mod- 



ulo a non-zero multiplicative term, one concludes from equation (|5.3|) that 

X 



u 



= c(t,u-)%£r, with c(t, in- 
putting this equation in equation (5.3), one finds 



a function depending on all variables except 



c(t, u. 



d 3 H d 3 H 



du 3 , ' didn't 



c(t, u. 



u+ 



<9|logP p 



d_ 

di 



d 2 E \og¥ v 



0. 



implying c(t, w_) = for all and thus H(t;u + ,u_) = \ogF v satisfies 
the equation X = 0, provided the Wronskian {<9f; logP p , -^d E \og¥' p } g ^ 0. 
This will be shown in the next section, using functional theoretical arguments. 
This ends the proof of Theorem 1.2 except for the "final condition", which 
will be shown in [TJ. ■ 



6 An estimate for the Wronskian 

Proposition 6.1 The Wronskian 



^llogP p ,9llogP p 



E 



is a non-zero function. 



In order to prove this proposition, we need a number of lemmas; the 
proofs will be functional-theoretical and rely on the techniques and on some 
of the formulae in [26]. Given the Pearcey kernel K v ', one defines, for a given 
set E = U^ =1 [a 2 fc_i, a,2k] C K, the following^ 



K V E := K V X E , I + R:=(I- KlY 1 =: p(x, y) (6.1) 
and thus one has 

(I-Kl)^Kl = Kl + (Elf + ... = R. (6.2) 



Also, from (1.2) and the differential equations (1.4) for p(x) and q(y), it 
follows that 

^ + -jj^X p (x, V )=p{z)q(y). (6.3) 



9 Given a kernel, viewed as an operator, the equality = refers to the corresponding 
kernel 
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Note that for a general kernel L, one has 

Define the functions 

V := (I - K^-'p, q:={I-K v E T )- l q (6.5) 

u := (I — K E )~ 1 p(a)q(a)da = (p(a), q(a)X E (a)). (6.6) 



E 



Lemma 6.2 For a disjoint union E = Ufc=i[ a 2fc-i? a 2k], one has the follow- 
ing identity~\ 



9llogP p = d|logdet(I - K V E ) = d E u = ^(-l) fc p(a fc )g(a fc ). 

k 

Proof: At first notice that 
[D,K%[ = [D,K?X B \ = (^ + |)(^) (6-T) 

= p(x)q(y)X E (y) - ^(-l) fc i^(x, a k )5(y - a k ), 



from which one deduces, using notations (|6.1|) and (]6.2|), 

<9x dy 



9 + = [D,{I-K^ 



- ^(-l) fc (/ - KlY'K^ix, a k )5(y - a k )(I - K^)' 1 
k 

p{x)q{y)X E {y) - ^(-l) fc i?(x, a k )p(a k , y). (6.8) 



w (f,9) ■= f M f(x)g(x)dx. 

11 Remembering Oe — Y^i afr- 
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and 



_d_ 
dai 



R{x,y) 



_d_ 
dai 



(I + R) 







{i -Kir 



da k 

(I - K v X E y x K p ^{I - K P X 



i v dX E tsp-y_\-i 



da k 

R{x,z){5{z-a k ){-l) k )p{z,y) 
(-l) k R(x,a k )p(a k ,y). 



(6.9) 



Then combining (|6.8|) and (|6.9|), one finds 

' d_ d_ 
dx dy 



k 



d 



[D,{i-Ki)- i \+Y.-^ R ^y) 



= P(x)q(y)X E (y), 



and hence, setting x — y — dj, the total derivative becomes 



(6.10) 



(6.11) 



One then computes the derivative of u, as defined in (]6.6|), with respect to 

ions j 

{I-K v E y x V ,qX E 



a k : (of course, the functions p and q do not involve the interval E) 

du d 



da k 



da k 







(—(I-K'P)- 1 )p,qX E \ + (p,q 



dX L 



'-l) k (R(x,a h )p(a k ),qX E ) + (p,qd(y - a k )(-l) k ), using (6.9) 



-l) fc (p(a k )(R(x,a k ),qX E ) + p{a k )q(a k )) 

-l) k p(a k )((I + R T )q(a k ) - q(a k ) + q(a k )) 
-l) k p(a k )q(a k ). 



(6.12) 
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Combining (6.11) and ( 6.12[ ) yields 



du 
daj 



d 



E _ ( _ imaj , aj . ; 



and then summing with respect to j 



12 



dai 



- £ 



da i 



logdet(/-7\~£ 



d% log det(J - Iff) = d\ log P 



which, together with (1.2), establishes Lemma 6.2 



Lemma 6.3 Given E = [x,y], the following estimates hold 
d d 



dx dy 



du 
di 



(y - x)(p(x)q(x))' + 0(y - x) 2 
-{y — x)(p"q — pq")(x) + 0(y — x) 



Proof: Using the fact that, for a small interval E, the integral J y has order 
x — y, using R(a, y) — R(a, x) = 0(y — x), one deduces from the formula of 



Lemma 6.2, (remember the definitions (6.5) of p and q) 
d E u = p(y)q(y) - p(x)q(x) 



y \ / rv 

P(y)+ I R(y,a)p{a)da) ( q(y) + / R(a,y)q(a)da 



y \ / rv 

p{x) + I R(x, a)p(a)da J I q{x) + / R{a 1 x)q{a)da 

2 



p(y)q(y) - p(x)q(x) + o(y - x) 

(y - x)(p(x)q(x))' + 0(y - x) 2 . 



12 Here one uses identity (1.1) in [55] 
d 



da 3 



logdet(J-Xl)- 1 = {-l) j+1 R{a,,a 3 ). 
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Using the heat equations (1.5) satisfied by p, q and the PDE (1.6) for K v , 
one checks 

2f = 2>-AS>-'»»> 



(/ - K 



V 



'at 



+((I-K^2% q ) + ((I-K^p,2 



dq 
di 



dy{-p'{x)q{y) +p{x)q'{y))p{y),({I - K E ) g) (ac 



-(/(*), ((/ - Kfy l q) (x)) + ((/ - Kl)- l p{x), q"(x] 



u(x)(-(p',q) + (q',p)) - (p",q) + (p,q"), using (p,q) = (p,q) = u, 

0(y-x) 2 -(p",q} + (p,q") 

— (y — x)(p"(x)q(x) — p(x)q"(x)) + 0(y — x) 2 , 



thus ending the proof of Lemma 6.3 



Proof of Proposition 6J_: One computes the following Wronskian and expand 
for small y — x, given E = (x,y). Indeed, from Oeu = d\ log (see Lemma 



6.2), the estimates of Lemma 6.3 and the the differential equations (1.4) for 



p and q, one computes 



d 



dt 



log 


p p ,9i 






d E 




ill 


-x) 2 



! 



d E 



({(p"q-pq"y,(pq)"} x + 0(y-x)) 



(y - x) 2 ( {-t(p"q - pq") + 3x(pq)' + 2pq)(pq)" + 0(y - x) 
*r{t(p'q - pq') - 2xpq + p"q' - p'q")(t(pq)' + 3(p'q')') 
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with the coefficient of (y — x) 2 /2, for x = t = being equal to 



2pq(pq)" — 3(p'q')'(p'q" — p"q') ^ 



which is nonzero, ending the proof of Proposition 6.1| 



7 Steepest descent analysis and replica dual- 
ity 



In this section, we emphasize the central role played by the spectral curve 
(or Pastur equation [21]) in the steepest descent analysis used to prove the 
universal behavior of the kernel as N — > oo. More precisely, we point out 
that the kernel used for the steepest descent analysis takes a very universal 
form in the different problems of non-intersecting Brownian motions studied 
so far. We further give show that the study of the spectral curve associated 
to the considered problem gives the universal limiting kernel obtained after 
a proper rescaling of the variable around a singularity of the problem: we 
give a "physical meaning" to the computations performed in the preceding 
sections as well as a way to generalize it to more complicated problems. 

In a first part, we show how this spectral curve arises in an integral 
representation of the kernel in the case of the matrix model in an external 
field. We then show how the expression of this kernel in terms of the spectral 
curves exhibits a universal behavior in the large matrix limit depending on 
the local properties of this curve. We finally apply this procedure to prove 
the appearance of the Pearcy and Airy kernels in the context described in 
the previous sections. 

7.1 From Hermitian matrix integrals to double con- 
tour integrals 

In this section, we derive the a double integral representation of the kernel by 
using the replica formulation introduced by Brezin and Hikami [H] . In order 
to make this paper self-contained, we show this duality by simply performing 
gaussian integrals. 

Let us consider the partition function 




(7.1) 
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where one integrates over hermitian matrices M of size N x N and A is a 



deterministic diagonal matrix, with arbitrary k of the form 



A := diag(ai, . . . , a{,'a 2 , . . . , a 2 \ ...,a k ,..., a k ). (7.2) 

Diagonalizing the matrix M and using the HCIZ formula [To! [17] , one is 
left with the integration over the eigenvalues (xi,X2, . . . ,xn) of M: 

f A(x) - N Y.(i- x ^) 

Z(A)= / n^XTY e 1 • ( 7 - 3 ) 

./V A(o) 

In order to compute the density of state -Ri(A) and the £-point correlation 
functions Rg defined by 



\i=l I 



R k (\ l ,\ 2 ,...,\z) = — ( J | Tr* A f .]- \f) ). 17. ii 

where the average is taken with respect to the probability measure 

1 « A (x) 

we consider their "Fourier" transforms 

U l (t lt U 1 ... i t l ) = (j[Tce^ M J (7.6) 

and, in particular, one gets the Fourier transform of the two points correlation 
function: 



N , / JV \ - N 



r 2 



ai,a 2 =l \i=l 



(7.7) 
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The previous section considers the particular case k = 2. 
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One can now integrate the variables Xj by noting that 



N 



N 



J \ Y[dxA A(x)e ^ 



2 / j J 

A(6)e j=1 



(71 



and expanding A(x) = Yli^j( x i ~ x j)'- 
U 2 {tiM) 



N 



AM it 1 a c , 1 +it2a ct2 ^2^-t 1 t 2 8 ai>c , 2 



l+*3 



X 



\<Km<N 



n ^ 



(7.9) 



l<l<m<N 

One can see that this can be written as a double contour integral 
U 2 (tuh) 

dudv 



e~ N 2 



tit 



1 L 2 



e 2 



(2z?r) 



3 ATi(t lU +t 2U ) (u-v + ih -it 2 )(u-v) 



X 



- LJ -\ u- a k J \ v - a k 



(u — v + it\)(u — v — it 2 ) 



it 2 



M 2 



dudv 



e Ni(tiu+t 2 v) / ]_ 

(2i7r) 2 \ (w — v + iti)(u — v — it 2 ) 



x 



x 



u-a k J\ v - a k ' 



it 2 



(7.10) 



where the integration contours encircle all the eigenvalues a k and the pole 
v = u — %t\ 



ii 



We can now go back to the correlation function 

dt\dt 2 



R 2 (X,ii) 



co poo 



CO J —CO 



47T 2 



N ^ x+t ^U(t 1 ,t 2 



(7.11) 



14 



see for instance [9] for more details around eq.(2-20) and eq.(4-40). 
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By first integrating on ti and t 2 with the shifts t\ — > t\ — iu and t 2 ^ t 2 — iu, 
this latter equation reads: 

R 2 (X, fi) = K N (X, X)K N (ji, fi) - K N (fi, X)K N (X, (7.12) 

where the kernel is defined by 

'MX,,) = / £ I £ ft ( —) -^e-"^^) (7.13) 
J 2vr / 2?vr JJL y v _ afc j v _ ^ 

where the integration contour for v goes around all the points a k and the 
integration for t is parallel to the real axis and avoids the v contour. More- 
over, it can be derived in a very similar way that any /c-point function can 
be written as the Fredholm determinant: 

R k (x u ...,x k ) = det [K N (xi, ^)]^ =1 • (7.14) 

By Wick rotating the integration variable t — > it, one gets 



K N (X,fi) = 




(7.15) 



where the integration contour for t is now parallel to the imaginary axis. One 
can then rewrite it under a more factorized form: 

K N (X, fx) =1^1 ^ e -"W*«>-W>> J_ (7.16) 
J 2m J 2m v — t 

where 

2 ^ 

S(x,y) = V — - xy + ^e;m(y - a t ) (7.17) 

with the "fraction numbers" given by 

e j: =ViV. (7.18) 

The first step in the steepest descent analysis of this kernel is to look for 
the stationary points of the exponent in the integrand, i.e. we look for y as 
a function of x solution of the equation 

k 

d y S(x, y)=y- X + J2 = (7.19) 

i=i y ai 
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which is nothing but the equation of the classical spectral curve introduced in 



the preceding sections (see 3.18) and in the general study of the one matrix 
model in an external field. 

Remark: If the external matrix A is highly degenerated, the ej are fixed and do 
not depend in N. Thus, the action S(x, y) does not depend on N in this case. 

Remark: The density of states is given by p(A) = Kn(\, A) and its derivative wrt 
A can be factorized: 

where 

0(A) = f !^ e ^(V) and m= I^L e -NS(x,v) (7 21) 

Jim 27r / 2«7r 

7.2 Saddle points, spectral curve and universality 

Let us now forget about the matrix model and consider a general kernel of 
the form 

K N (x,fj) = (-if 1^1 £L e -ms M - S{ x, m i (722) 

J Z7l J 2,171 V — t 

where the derivative of the action S(x,y): 

is a rational function which can be written as the ratio of two polynomials 
in both variables. Then the locus of stationary points of S is given by an 
algebraic equation 

£(x,y) = 0, (7.24) 

referred to as the spectral curve in the sequel. Let us study some of its 
properties necessary to classify the different universal behavior of Kn(\,h) 
as N — > oo. 

Generically, the equation £(x,y) = has d y distinct solutions in y for a 
given value of x. Let us denote them Yi(x) as functions of x: 

dy 

S(x,y) = g dy (x)H(y-Y i (x)) (7.25) 

i=l 
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where gd y (x) is the leading coefficient of £(x,y) as a polynomial in y. How- 
ever, there exists finitely many branch points Xi such that £{x,y) = has 
a double zero, i.e. two solutions Yj(xi) = Yi(xi) coincide. One can also 
characterize them by the property: 

d y £(xi,y)\ y :=Y j{xi) = 0- (7-26) 

More generally, a £th order branch point x!p is defined by 

97 £(4\v) Y . = 0, for m < £, with ^f, y) ^ 0. 

(7.27) 

Other types of singularities might occur, but we will not refer to them in this 
paper. 

Remark: A (£ + l)th order branch point can be obtained when a £th order branch 
point and a simple branch point merge. 

In the next sections, we show that, as iV — > oo, the kernel K N (X,/i) 
has a universal behavior for A and /i approaching the same point with an 
appropriate scaling (depending on iV and the singular behavior of the spectral 
curve at this point). We study the first singularities and show this universal 
behavior by a local analysis of the spectral curve whereas the usual Riemann- 
Hilbert study [?] involves a global analysis of the latter. 



7.3 Steepest descent analysis 

The asymptotic of this kernel when N — > oo exhibits different regimes de- 
pending on the value of its argument. More precisely, we get a universal 
kernel associated to the neighborhood of any point Xq of the spectral curve: 
if Xq is a generic point of the spectral curve, we get the sine kernel; if Xq is 
a simple branch point, we get the Airy kernel; if xq is a higher order sin- 
gularity of the curve, we get a universal kernel associated to this particular 
singularity. 

In any case, we proceed with the same steps: we choose a point Xq. We 
then focus on its neighborhood by a change of variable consistent with the 
singularity of the curve at xq. We expand the exponent in the kernel around 
this point by simply writing a Taylor series expansion. 
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7.3.1 Simple branch point: the Airy kernel 



Consider a one-parameter family of algebraic functions S(x, y\t), parametrized 
by t, such that there exists a critical point (x c , y c , t c ) satisfying 



Sy(x c . y c i c j Syy{x c: y c \t c ) 0, 

) 7^ 0, S yt (x c ,y c \t c ) ^ 



(7.28) 
(7.29) 

It means that the spectral curve S(x,y\t c ) has a simple branch point at 
(x c , y c ). One shows that in the neighborhood of this branch point, the kernel 
can be rescaled in such a way that it converges to the Airy kernel as N — > oo. 
Consider the changes of variables allowing to focus on the neighborhood of 
the critical point: 



t 


= t c - 








x 




+ 






y 


:= Vc 


_|_ a y Y 


NT. 



I3 X X 

2 

AT3 



(7.30) 



and expand S(x,y\t) around this critical point using (7.30) as N — > oo. One 
then expands S(x,y\t) and S(x,y\t) in aTaylor series in t, x and y, thus 
obtaining for the expression in the integrand of the kernel: 

N (S(x,y;t) - S(x,y;t)) 

= iVs \ a y [a x S xy (x c , y c ; t c ) + a t TS ty (x c , y c ; t c )) (Y - Y) + (3 x XS x (x c , y c ; t c ] 

3 ~3 a y 
~\~(Y Y ) Syyy(x c , y c ,tc) 



6 



+ {Y 2 - Y 2 



a 



p [a x S xyy (x c , y c ; t c ) + a t TS ty y(x Cj y c ; t c )\ 



+ (XY - XY)f3 x a y S xy (x c , y c ; t, 
+(X-X)(3 X 

+ (Y - Y)Oy 

+0(iV-s); 
which we write for brevity as 
N(S(x,y;t)-S(x,y;t)) = 



a x S xx (x c , y c ; t c ) + a t TS tx (x c , y c ; t c )] 

S xxy (x c , y c \ t c ) + ol x ol{F St X y(x c , y c \ t c ) H — T Stt y (x c , y c \ t c ) 



(7-31) 



AT3 



ai (y — y I +a 2 X-X 



3 I i / v2 y 2 



+« 5 xy - xy 



« 6 ( y - y 
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+ « 3 (y 3 - y 5 ) + n, (V 

h a 7 ( X - X 



(7.32) 



for coefficients {oj}J =1 functions of the scaling parameters a t , a x , (3 X and a y . 

The coefficients 02 and a.j can be eliminated by conjugation of the kernel 
leaving the Fredholm determinants invariant. In order to recover the Airy 
kernel, one has to fix the remaining coefficients by 

oti = , a3 = I/3 , 05 = — I and 0:6 — 0:4 = 0. 

(7.33) 

a y and f3 x are determined by the constraints on a 3 and o 5 respectively, while 
the first equation gives and the last one determines a 
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7.3.2 Double branch point: the Pearcey kernel 

Let us now consider an algebraic a function S(x,y\t) such that there exists 
a critical point (x c , y c , t c ) satisfying 

SyijEci Vc\tc) Syy{X C: y c \tc) S yyy (^X C , y C \t C ) 0, (7.34) 

Syyyy(x c ,y c \t c ) ^ 0, Sy X {x c ^y c \td) 7^ 0, Syt(x c , y c \t c ) 0. (7.35) 

It means that the spectral curve S(x,y\t c ) has a double branch point at 
(x c , y c ). In the neighborhood of this critical point, one can rescale the kernel 
in such a way that it converges to the Pearcey kernel as iV — > oo. To this 
effect, consider the changes of variables in the neighborhood of the critical 
point: 

x := x c + + ^# (7.36) 

AT2 ATI V ) 

V'=Vc + -*r 

ATI 



and expand S(x,y\t) around the critical point using (7.36) as N — > oo. 

N (S(x,y;t) - S(5t,y;t)) = 



: ATI 

at 



(3 x S x (x c ,y c ;t c )(X - X) + a y (Y - Y) (a x S xy (x c ,y c ;t c ) + a t TS ty (x c ,y c ;t c )) 



24 ^ yyyy ^ Cc ' 1 ^ c )(^ ) 
a 2 

+ Y [a t TS t yy(x c , y c ; t c ) + a x S xyy (x c , y c ; t c )] (Y 2 - Y 2 ) 

+(3 x a y S xy (x c ,y c ;t c )(YX - YX) + 0(N~*). (7.37) 

Remark: One can see that (as in the previous case) all the terms in this expression, 
except the scaling factor f3 x S x (x c , y c ; t c )(X — X), depend only on derivatives of S 
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wrt y. This shows that the kernel depends only on the spectral curve £(x, y\t) := 
S y (x,y\t). 

One can see that the variables to be integrated Y and Y appear only in 
terms which do not blow up as iV — > oo except one which is proportional to 
Ni. One can get rid of this term by fine-tuning the coefficients of the change 
of variable. Indeed, imposing the constraint on 

a x S xy (x c , y c ; t c ) + a t TS ty (x c , y c ; t c ) = (7.38) 

eliminates this term. One can finally normalize the scaling coefficients a y , 
a x and (3 X respectively by setting 

a t 1 

2^^yyyy^ c ^ PxayS X y(^x c ,y c ,t c ) 1 (7.39) 



2 

OL X OLy 



a t T 



Sfyy{Xci Va tcj ~\~ S X yy{Xci Vcj t c 



T 

2' 



(7.40) 



a x 

yielding the quartic exponent appearing in the Pearcey kernel: 

N(S(x,y;t)-S(x,y;t)) 
= N l if3 x S x (x c ,y c ;t c )(X -X) 

+ I(y4 _ _ L(y 2 - Y 2 ) + (YX - YX) + O(Art), (7.41) 

where the term in (X — X) can be eliminated by conjugation of the kernel. 

7.3.3 kth order branch point 

More generally, consider a point (x c ,y c ;t c ) of the spectral curve £(x c ,y c ;t c ) 
where^l 

cT£(x c , y c ; t c ) = 0, for all m < I, with d y +1 S(x c , y c ; t c ) ^ 0, (7.42) 

while still assuming that dtS and d x 8 do not vanish at this critical point. 

Let us use this example to explain how one can guess the rescaling to 
obtain the desired universal kernel. Generically, one looks for a rescaling of 
the form: 

+ ._+_)_ atT 

x:=x c + fi (7.43) 

, a y Y 

y • yc ~t~ N i y 



15 Remark that the Airy case corresponds to I = 1 and the Pearcey case to I = 2. 
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for some critical exponents 7$, 7 X and 7^ to be determined. One then writes 
the Taylor expansion of the action S around the critical point using this 
rescaling: this expansion can be seen as a series in iV -7 for some exponent 7 
depending on the rescaling. Since the action is multiplied by N in the kernel 
and integrating over the variable y, one needs that all the terms depending 
on Y in this expansion are of order at most as N — > 00. In particular, 
considering the Taylor series with respect to y: 

S(x c ,y;t c ) = S(x c , y c ; t c )+- ■ ■+ d l y +1 S(x c , y c ; t c )N^ l+2 ^+o (jV~^) , 

one must impose (I + 2)^ y = 1 for this contribution to appear, thus fixing 
the critical exponent for y: 

7. = (7.44) 

On the other hand, one expects x and y to couple in the exponent of the 
limiting integrand (otherwise, it would only give a simple multiplicative fac- 
tor). Thus, the first contribution of a mixed derivative wrt to x and y of the 
action must be of order which imposes: 

7, = 1-7, = ^- (7.45) 

One can remark that we could stop the procedure here and not rescale 
the time which could remain decoupled from x and y. But one could also 
try to couple it with y by taking a t 7^ 0. If one wants this coupling to be 
different from the one between x and y, one must impose 74 < 7 X and one 
gets, in the simplest case 



16 



7t = l-27 y = 77^. (7.46) 

By doing so, one has coupled t to y. But one has also introduced a divergent 
term coming from the coefficient of Sty One can compensate this term by 
completing the change of variable for x and writing: 



* :=t c + — r 

iyi+2 Ni+2 y ' 

y :=y c + ~^t- 

Nl+2 



16 This means that one wants to have a coupling TY 2 , but one could also obtain in a 
more complicated way higher order times by looking for couplings of the form TY k with 
k > 1. 
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Doing so, one finally gets the associated kernel in terms of the rescaled 
variables: 

K k (X,X>) = jdY J dY'e^- T -t^ Y e-^- T -^ +x ' Y ' y^y,. (7.48) 

Remark: One omits a detailed proof here, although the details can be filled in. 
Indeed, in order to do so, one carefully studies the integration contours and the 
normalization, which implies a rescaling of the kernel as well as conjugation. It 
is done in the preceding sections for the Airy and Pearcey cases for two ending 
points. 

7.4 General case 

In a more general case, i.e. the two arguments x and y of K N approaching 
any singular or non-singular point of the algebraic curve, one can use the 
same study leading to a different universal limiting kernel associated to each 
type of singularity. Let us summarize the procedure one has to follow to 
obtain the universal kernel associated to a given singularity: 

• First compute the multiple derivatives of the action S(x, y; t) with re- 
spect to x, y and t at the considered critical point and fix which are 
the first non-vanishing derivatives; 

• Fix the critical exponent in the rescaling by studying the large N be- 
havior of the Taylor expansion of the action around the critical point, 
under consideration; 

• Normalize the change of variables and compute the rescaling of the 
kernel through the Taylor expansion of the action of the action around 
the critical point: one then gets the universal limiting kernel; 

• Finally, carefully study the integration contours to check that they give 
a right path for the steepest descent analysis. 

7.5 Example: back to the matrix model and non-intersecting 
Brownian motions 

Let us now apply the preceding procedure to the case of the random matrix 
model in an external field studied in the preceding sections: one considers N 
Brownian motions starting from and going to two ending points at t — 1. 
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As reminded in section[2|of this paper, the kernel, given in Eq. (2.2), can 



idu f dv 



easily be taken to the form of Eq. (|7.22|) 

H n (fl,X;t) = 
by using the rescaling 



■n(S(p,,v;t)-S(X,u;t)) _ 



2vr / 2in 



v — u 



vt 



ut 



at 



u 



C\/n 



C\/n 



b :-- 



bt 



C\/n 



C\/n 



/i :-- 



x 



and A := 



n 



with 



and the action 



(7.49) 



y_ 

fn 

(7.50) 



S(x, u; t) :- 



u 



u 



x — h p \og(u — a) + (1 — p) \og(u — b) 

c 



(7.51) 



which is nothing but the function S(u,x;t) = F(u) of Eq. (3.18) studied in 
section |3| The study of this section, and particularly Eq. (3.20), states that 



(7.52) 



S u (xo, uq\ to) — S uu (xo, uo; to) — S uuu (xo, wo; to) — 
at the critical branch point (xq,uo) of the spectral curve 



S u (x,u;t ) = u 



x 



p 1 
— + 



P 



+ 

c u-a u-b 







(7.53) 



where one notes c = c(t ), a = ^ and 60 = ^- Moreover, using the 
notations and computations of section [3j one has: 



'S'ua^O; Uq] to) 



Co 



(x ,u ;t ) 
4! 



q+ 1 



Aq 



(7.54) 



and 

S u t{xoi uq; to 



Wl - 2t c 



+ 



pa cl 



+ 



;i -p)boc 2 o 



2c t (l - t ) tl{\ - t ) 2 {u - So) 2 - t ) 2 (u - b o y ' 

(7.55) 

This implies we are in the case of a double branch point studied in section 



7.3.2 and thus, with the right rescaling, the kernel converges to the Pearcey 



kernel as n — > 00. Let us now check the conditions to obtain precisely 
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the Pearcey kernel. From the study of section 7.3.2, one must consider the 
rescaling 

' • ■ to 



t :=t + ^ 



x := x + + ^# 
u := u + 

with the conditions 

Oi x Sxu(xo, Wo! to) + OltTSui(x , u ; t ) = 

<*u 1 

^'<S'uuuu(^05 M 0i to) = 



(7.56) 



(7.57) 
(7.58) 
(7.59) 

f3 x a u Sux(xo,u ;to) = -1. (7.60) 

Let us check that these conditions agree with the rescalling Eq. ( 3.17[ ). 
From conditions (7.58) and (7.60), one gets (/x was defined in Eq. (1.10)) 



and 



Of T 

-y [^^^(xo, «o; to) + O^S^i^O, «o; to)] = g 



a u — — , /3 X — Co/i, 



(7.61) 



whereas Eq. (7.59) and Eq. (7.57) give 



;i - 1 ) 2 ^ 



pap 



+ 



(l-p)feo 



(u -io) 3 (mo-6o) 3 



g 2 - g + 1 



2c 2 /! 2 



(7.62) 



and Eq. (|7.57|) states that 

x (l-to)(l-2t ) 



a, 



2c n 



+ 



pa cl 



[l-p)~b c 2 
t Q (u - a ) 2 ' t (u - b ) 2 



+ 



pap 



+ 



(l-p)feo 



(n -a ) 3 1 (« -M 3 

It is then easily checked by plugging in the values of x , t and u in terms 
of q (with Maple for example) that this rescaling coincides with Eq. (3.17), 
that is to say, the rescaling considered in the preceding part. 
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7.6 Application: Brownian bridges from one point to 
k points 

We can also use this analysis to study more general statistical systems. Let 
us now consider N non-intersecting Brownian motions stating from at time 
t = and ending at k points dj by groups of rii particles where Hi = ejiV. It is 
a classical result that their probability measure at a given time t is given by 
the probability measure of the eigenvalues of an Hermitian random matrix 
in an external field with external matrix A(t) whose eigenvalues are given by 

Then, the spectral curve is given by: 

£(,,„)=„-, + ^-5^ = (7.63) 



i=l 



and has genus zero. One can straightforwardly (just solving the equation in 
x) find a rational parametrization under the form: 

x{x ) = l + ± J -± m , y {z) = z (7.64) 

and the branch points are given by 

^ ) = 1 -g(73^ = °- (? ' 65) 

For a generic time t, one has 2k branch points, 21 of them {z2i-i, Z2i)\ =l 
being real and lying in the so-called physical sheet of the spectral curve 
(see [HI [8j [23] for an extensive study of the spectral curve). Thus, by he 
study of section [3] and |7.3[ one can conclude that the kernel converges to the 



Airy kernel in the neighborhood of these real branch points Zj recovering the 
results of [2"3] . 

Now, as the time decreases from 1 to 0, some real branch points Z2i(t) and 
Z2i+i{t) merge for some critical value of the time t c . In the Brownian motion 
setting, it correspond to the times when one big group of particles splits into 
two smaller. In terms of the family of spectral curves parameterized by the 
time t, it corresponds to the merging of two cuts. In the neighborhood of 
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this double branch point, thanks to the study of section 7.3, one can rescale 



the kernel in such a way that it converges to the Pearcey kernel as N — > oo. 



This generalizes the result of th. 1.1 to any cusp in this kind of processes 



and recovers the former results of B3 



It is then natural to ask the question: Can we have higher order singu- 
larity in these processes? When the e^'s are independent of iV and the a^'s 
real, the answer is no. 

Indeed, it amounts to knowing what is the highest order possible for a 



real root of Eq. (7.65 ). This problem can be rephrased as knowing the higher 



order possible for a real root of the equation 



(7.66) 



2t 
1-t 



with the constraint Y2i e « = 1 an d the rescaled time evolution T 
[0,oo]. ^ 

One can prove that the real roots of this equation are at most double. 
For this purpose, let us follow the evlution of the roots as t decreases from 1 
to 0, i.e. for T going from oo to 0. 

For T large, this equation has obviously 2k simple real roots Z{ located 
around their T — ► oo value, i.e. \z 2 %-\ — a%\ C 1 and \z 2 i — a<| <C 1 with 
Zn-\ < cii < z 2i . Now, for any real solution of this equation, one can compute 

This derivative does not change sign as long as Zj does not cross any eij. It 
means that z 2 i (resp. z 2i+ i) keeps on going from to a i+1 (resp. from a i+ \ to 
dj) as T decreases, unless it reaches another real root. Following this process 
one sees that the two real roots z 2 i and z 2 i + \ meet for some critical time T c 
giving birth to a double real root of this equation. In this first part of the 
process, one can thus only encounter double real roots. 

Let us keep on decreasing time. For T < T c and close to it, the double 
root gives rise to two simple complex conjugated roots. Let us thus now 
consider a simple complex root z = r + i6. For 9 close to 0, one can compute 

de i . 
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to first order in 9. The complex roots are thus repelled by the real axis and 
cannot thus give birth to real roots. In this second part of the process, one do 
not have real roots anymore. The only multiple real roots are thus obtained 
when a,i < z 2i z 2i+1 < a i+1 . 
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